PLh ' Table 1. OB Stars 



(N ; star SpType Ref 

HD 157955 B9 IV 1 

> HD 169033 B8 IV-Ve 1 

HD 183143 B7 la 2 

BD+24 3866 08.5 lb 3 

^ HD 190429 04 If+ 4 

O , HD 190603 B1.5 Ia+ 5 

I^ ■ HD 190864 06.5 HI(f) 6 

^ HD 198183 B5/6 V 1 

^ . HD 198478 B2.5 la 5 

C^: HD 199579 06 V((f)) 4 

Q HD 201254 B3 V 1 

^ ■ '^ SpType references: 1.) Bright Star Catalog, Hoffleit (1982) 2.) Morgan, Code, & Whitford (1955); 3.) Conti 






a 



Alschuler (1972); 4.) Walborn (1973); 5.) Lennon, Dufton, & Fitzsimmons (1992); 6.) Walborn (1972); 

Notes to Table 1. 



X 

^ ■ All spectra have a two pixel resolution of A/AA ~ 570, except HD157955 and HD169033 which have A/AA « 525; see 
<^ ' text. 
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Table 2. Hydrogen and Helium Equivalent Widths (A) 





Line 


A i^m) 


HD190429 


HD199579 


HD19G864 


BD+24 3866 


HD190603 


HD198478 


HD2G1254 


HD198183 


H 


Pa7 


1.0938 


2.0± 0.2 


3.9± 0.5 








4.7± 0.2 




7.1± 0.2 


r 


i-H . 


Pa/? 


1.2818 




5.6± 0.7 


4.0± 1.0 


4.4± 1.0 


4.8± 0.7 


6.0± 0.8 


8.5± 0.8 


9.1±0.8 


C 


2237\ 


Brio 
He I 


1.7362 

1.7002 


1.4± 0.5 


1.8±0.5 


2.4± 0.8 


3.6± 0.4 


2.7± 0.9 
3.0± 0.4 


2.9± 1.0 
2.4± 0.4 


6.2± 0.8 
1.8±0.5 


7.3± 1.1 
1.6±0.6 


c 

\ 


Brll 
Brl2 


1.6807 
1.6407 












3.4± 0.9 
3.4± 0.9 


3.3± 1.0 
3.3± 1.0 


6.3± 0.8 
5.0± 0.8 


7.7± 1.1 
6.3± 1.1 


c 


o, 


Brl3 
Brl4 
Bil5 


1.6109 
1.5880 
1.5701 












3.0± 0.9 
1.3± 0.9 
1.8± 0.9 


2.7± 1.0 
1.6± 1.0 
2.2± 1.0 


4.1± 0.8 
2.2± 0.8 
2.1±0.8 


4.5± 1.1 
2.4± 1.1 
2.1± 1.1 


c 




Brl6 
Brl7 


1.5557 
1.5439 












1.1± 0.9 
0.4± 0.9 


1.2± 1.0 
1.3± 1.0 


0.5± 0.8 


0.6± 1.1 




o 


Brl8 


1.5342 












0.9± 0.9 


1.0± 1.0 








• • 


Brl9 
Br20 


1.5261 
1.5192 












0.8± 0.9 


0.8± 1.0 








Br21 


1.5133 
















... 







Notes to Table 2. 

Hydrogen and Helium absorption strengths are in A. The line wavelengths are given in air (see text). Lines were measured by fittin 
are the formal one sigma errors based on the rms scatter in the continuum and do not reflect the line fit. 



0^ 



PLh 



en 
(N 



cd 



i7— Band Spectroscopic Classification of OB Stars 

R. D. Blum\ T. M. Ramond, R S. Conti 

JILA, University of Colorado 

Campus Box 440, Boulder, CO, 80309 

rblum@ casa . color ado . edu 

ramond@casa.colorado.edu 

pconti@jila.colorado.edu 



-^ ■ D. F. Figer 

O ' Division of Astronomy, Department of Physics & Astronomy, 

University of California, Los Angeles, CA, 90095 
^ ■ figer@astro.ucla.edu 



. , K. Sellgren 

C^ I Department of Astronomy, The Ohio State University 



174 W. 18th Ave., Columbus, OH, 43210 



CN I sellgren@payne.mps.ohio-state.edu 

t^^ ! accepted for publication in the AJ 

^- ABSTRACT 

O 

H 

^ \ We present a new spectroscopic classification for OB stars based on if —band 

(1.5 /im to 1.8 /im) observations of a sample of stars with optical spectral types. 

Our initial sample of nine stars demonstrates that the combination of He I 

rS ' 1.7002 fiui and H Brackett series absorption can be used to determine spectral 

c3 I types for stars between ~ 04 and B7 (to within ~ ± 2 sub-types). We find 

that the Brackett series exhibits luminosity effects similar to the Balmer series 

for the B stars. This classification scheme will be useful in studies of optically 

obscured high mass star forming regions. In addition, we present spectra for 

the OB stars near 1.1 /im and 1.3 ^m which may be of use in analyzing their 

atmospheres and winds. 



Subject headings: infrared: stars — stars: early-type — stars: fundamental 
parameters 
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1. INTRODUCTION 



OB stars are massive and thus short hved. Because they have short hves, they will 
be confined to regions relatively close to their birthplaces and will be found close to 
the Galactic plane. The youngest massive stars will also still be in or near their dusty 
star forming environment. Due to the strong extinction from interstellar dust at optical 
wavelengths, OB stars at large distances from the sun, or in young star forming regions, 
are more easily studied at near infrared (1 — 5 fim) or longer wavelengths. Note that the 
extinction at V is about 10 times greater in magnitudes than at 2.2 /im, and for ^ 3.5 A 
^ 0.9 /im, the extinction can be represented by a power-law dependence: Ax / ^2.2^m = 
(A/2.2)-i-^ (Mathis 1990). 

The advent of sensitive infrared array detectors has made near infrared classification 
schemes commonplace for a wide range of stellar spectral types. For example, Kleinmann 
& Hall (1986), Greene & Meyer (1995), Ah et al. (1995), and Ramirez et al. (1997) have 
presented i^— band spectra of late-type stars useful for classification at 2 /xm. Langon 
& Rocca-Volmerange (1992) have presented H and K spectra of optically classified stars 
ranging in spectral type from 06 to M7 (only three OB stars were included). Eenens et al. 
(1991), Blum et al. (1995a), Tamblyn et al. (1996), Morris et al. (1996), and Figer et al. 
(1996, 1997 ) have presented ii'— band spectra of known Wolf-Rayet and other emission-line 
stars. A detailed classification scheme in the i^— band for early-type (OB) stars has been 
presented by Hanson & Conti (1994) and Hanson et al. (1996, hereafter HCR96). All these 
infrared classification spectra are now routinely applied to the study of obscured stars, 
typically in star forming regions and/or along sight lines with large interstellar obscuration 
by dust (e.g. Greene & Meyer 1995; Blum et al. 1995a, 6; Tamblyn et al. 1996; Figer et al. 
1996; Ah 1996; Blum et al. 1996). 

Classification of hot stars at near infrared wavelengths is challenging because the 
photospheric absorption features used in the classification (primarily H, He I, and He II) 
are generally weaker than the molecular and atomic (metal) lines used for late-type 
stars. In addition, there can also be strong circumstellar emission which dilutes the 
photospheric lines. It is this latter problem that motivates the present work. Since we 
expect circumstellar dust emission to dominate the iT— band spectra of some stars in star 
forming regions (e.g., in M17, Hanson & Conti 1995), a more robust classification scheme 
for high mass star forming regions will also utilize shorter wavelengths which may be less 
affected by circumstellar excess emission. Conversely, we wish to observe at the longest 
possible wavelengths shortward of K to overcome the extinction of intervening dust. In this 
paper, we present a preliminary classification of OB stars based on spectra between 1.5 /im 
and 1.8 /xm (i.e. in the if —band). We also present spectra of the same stars between 1 /xm 
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and 1.3 /iin which may prove useful in the analysis of massive stars atmospheres. 



2. OBSERVATIONS and DATA REDUCTION 

The OB stars were observed on the nights of 9, 10, and 13 June 1996 using the 
Ohio State Infrared Imager and Spectrometer (OSIRIS) on the Perkins 1.8-m telescope of 
the Ohio Wesleyan and Ohio State Universities at the Lowell Observatory. The Perkins 
telescope is located on Anderson Mesa near Flagstaff, Arizona. OSIRIS is more fully 
described by DePoy et al. (1993). The target stars were taken from the catalog of IICR96 
and are listed, along with their optical spectral types, in Table |I]. In addition to the target 
OB stars, we observed A-type and G-type stars for use in canceling telluric absorption 
features. We will refer to these stars as "atmospheric standards." 

OSIRIS was used in cross-dispersed mode which gives A/AA ~ 570 (2 pixels) while 
covering the J (Ao ~ 1.25 /xm, AA ~ 0.30 //m), H (Ao ~ 1.65 /im, AA ~ 0.37 yum), and 
K (Ao ~ 2.20 /im, AA ~ 0.50 /im) bands simultaneously. One star (HD183143) was also 
observed in long slit mode at J and H. Long slit mode is similar to the cross-dispersed 
mode (same spectral resolution and spatial scale) except only one wavelength band is 
acquired at a time and the long slit fills the entire spatial dimension on the array. Similar 
results were obtained from both sets of spectra; however, the long slit mode spectra are 
slightly higher S/N and are therefore presented here. In addition, five of the stars observed 
in cross-dispersed mode were also observed at A ~ 1.10 yum (AA ~ 0.25 /im), a bandpass 
which we call /', in long slit mode. Analysis of iJ— band OH sky lines results in linear 
dispersions of 9.70, 11.60, 14.53, and 19.37 A pix^^ at /', J, H, and K, respectively. The sky 
lines were identified using the list of Oliva & Origlia (1992). /', J, and K band dispersions 
follow from the ilf— band dispersion and the appropriate order number, m (m = 6, 5, 4, 3 
for /', J, if, K). In cross-dispersed mode, OSIRIS has a ~ 60" x 5" slit (1.5" pix~^). The 
target stars were observed at ^ 12 uniformly spaced positions along the slit. Combining 
spectra over a uniform grid along the slit greatly reduces systematic errors which may be 
introduced, for example, by scattered light in the dome flats and fringing (see Blum et al. 
1995a for a discussion of these problems in OSIRIS). The seeing varied throughout each 
night from 2" to 4". None of the nights were photometric; we observed through thin clouds. 
This will not affect our results since we are interested in the relative intensities and line 
strengths, not absolute fluxes. 
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All basic data reduction was accomplished using IRAF[|. Each image was flat-fielded 
using dome fiats and then sky subtracted using another image from the grid with the 
star displaced by several positions along the slit (cross-dispersed mode) or with a median 
combination sky image (long slit mode). Individual spectra were extracted from each 
program star image and atmospheric standard image using IRAF "APEXTRACT." 
Synthesized apertures ± 3 pixels wide were used. The entire grid of 1-d spectra for each 
star was then combined (after scaling). In the case of the long slit spectra, the individual 
1-d spectra were first shifted (0 to ~ 4 pixels) to account for anamorphic demagnification 
along the slit spatial dimension. 

The final spectra were obtained by ratioing the program stars with a atmospheric 
standard star which had first been corrected for intrinsic H absorption lines (P/? 1.281 /im, 
and Brackett series 1.51 /im to 2.165 fim). The atmospheric standards were observed in 
pairs of A-type and G-type stars. The H lines in the A stars were corrected in the following 
way. The majority (seven) of our program stars were corrected using an AO V and a G2 
V. In this case, the 02 star (BS 6847) was first corrected for H line absorption using the 
NOAO solar atlas (Livingston & Wallace 1991). Next a ratio of the A star (BS 7734) to 
the corrected G star was made. Since these two were taken at nearly the same airmass, the 
resulting ratio contains essentially only the H line spectrum of the A star and perhaps some 
"emission" lines due to metal lines in the G star. The H lines in this ratio were then fit with 
Gaussian profiles. The resultant line fits were used to correct the H absorption in the A 
star. The two remaining stars (HD183143 and BD+24 3866) were corrected using an A3 V 
(BS 7958) and a G9 III (BS 7760). We have no intrinsic spectrum to correct BS 7760 with, 
so the J and i^— band lines were corrected by eye. No correction was made for the if —band 
Brackett lines in BS 7760. Since BS 7760 was not corrected using a matching intrinsic H 
spectrum, we are less confident in the precise H line absorption resulting in HD183143 
and BD+24 3866; however, none of our primary results is sensitive to the accuracy of our 
measurement of H absorption in the program stars but only the appearance or absence of 
a definitive Brackett series. No telluric correction was made for the J' spectra as a cursory 
inspection of the spectra shows no strong telluric lines. We did obtain A and G star spectra 
for use in making a similar telluric correction as described above. Making this correction 
shows no difference in the final /' spectra of the OB stars at the few percent level and serves 
only to introduce more noise. 

We have also included two if —band spectra (HD 157955 and HD 169033) which were 



^IRAF is distributed by the National Optical Astronomy Observatories. 



-5- 



obtained as part of a program concentrating on i^— band spectra of stars in an inner Galaxy 
star forming region called the Quintuplet (Figer et al. 1996). These spectra have similar 
spectral (16.75 A pix~^) resolution as the OSIRIS data and spatial scale somewhat smaller 
(0.7" pix~^); see Figer et al. (1995). The telluric correction was made by ratioing the stars 
to a Quintuplet star (q3 in Moneti et al. 1994) and a dusty Wolf-Rayet star (WR118), 
neither of which have spectral lines near the He I or Brackett lines. These stars were 
observed at higher airmass (~ 2) than the other stars presented here (~ 1) and may not be 
corrected as well (the higher Brackett lines appear to suffer from incomplete correction). 



3. RESULTS and DISCUSSION 

3.1. The iJ-band Spectra 

The final if —band spectra are shown in Figure |I|. To highlight the absorption features, 
we have divided each spectrum by a low order fit to the continuum. We do not present our 
/C— band spectra here since spectra for these stars have been presented by IICR96 (at higher 
spectral resolution). We note that our lower resolution K—hand spectra reproduce the 
basic features of the HCR96 spectra. Line identifications used throughout this paper were 
taken from Weise et al. (1966, H, He I), and Garcia & Mack (1965, He II). All wavelengths 
used and/or quoted are in air. Line equivalent widths (Table |^) were measured by fitting 
Gaussian profiles using the LINER program (Pogge 1997). Uncertainties are the formal one 
sigma errors derived from the rms scatter in the continuum corresponding to a bandpass 
equal to the full width at zero intensity of the lines. The errors do not reflect the line fit. 



3.1.1. The Brackett Series 

The H lines are sensitive to the excitation in the stellar atmosphere. For example, the 
strength of the well-known Balmer lines increases with T^^ as more atoms are excited to 
the n=2 level. The Balmer lines reach peak strength at ~ 9000 K (Gray 1992) as Hydrogen 
becomes appreciably ionized at higher temperatures. The Brackett series excitation is 
similar to that for the Balmer series (10.8 eV compared to 10.2 eV), so we expect the lines 
to behave in a similar manner. This basic behavior with Tgfj (spectral type) is seen in 
Figure |l|. 

Inspection of Figure ^ suggests the Brackett lines behave similarly to the Balmer lines 
in early-type stars as first noted by Adams & Joy (1922, 1923): the width of the H lines 
decreases with stellar luminosity for the B stars; further, the last line visible in the series 
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increases with luminosity. Brl9 is confidently detected in the B supergiants while Brl5 or 
Brl6 is the last of the series for the B dwarfs. The measured line widths (FWHM ~ 60 A 
for Brll-Brl5) for the two B dwarf stars are roughly two times the instrumental resolution 
(FWHM = 29 A). This can be compared to a maximum of 44 A in the B supergiants for the 
same lines. These "luminosity" effects are well understood as resulting from the linear Stark 
effect (Hulbert 1924; Struve 1929; Gray 1992) which is sensitive to the electron pressure. 

In Figure |I] it can be noted that the upper Brackett lines are missing in the earliest 
0-type stars, although Br7 is present in these same stars (HCR96). There are, as yet, 
no published model predictions for the upper Brackett series lines, but we suspect their 
absence in the earliest 0-type stars, while Br7 is still present, is analogous to the case for 
the Balmer lines. LTE model predictions of Auer & Mihalas (1972) clearly underestimate 
H7 absorption compared to observations, while their non-LTE predictions provide a much 
better fit (Conti 1973). Peterson & Scholz (1971) demonstrate that the difference between 
the observations and LTE predictions is much less for the higher series member H8 than for 
H7. The upper series members of the Brackett lines probably are still close to LTE (similar 
to H8) and thus, quite weak due to the advanced ionization state of hydrogen. The lower 
series members, such as Br7, exhibit non-LTE behavior (similar to H7) and are thus still 
present. 



3.1.2. Hel and Hell 

A relatively strong line of He I (4-3) is visible at 1.7002 /im. To be useful for 
classification purposes, the wavelength range of interest must contain various lines sensitive 
to different excitation in the stellar atmosphere. In the present case, the Brackett series 
lines (as discussed above) and He I 1.7002 /im line are sensitive to lower and higher T^s, 
respectively, in the atmospheres of OB stars. In Figure |^, we plot the equivalent width of 
the He I (4-3) 1.7002 /xm line versus the spectral type (Table |l|). This line is at maximum 
strength in the late 0/early B types. In the earliest O stars, the Brackett series will be weak 
or absent (Figure |I]), whereas in the late B stars, the Brackett lines are stronger. These 
features, therefore, provide a rough spectral class in the if —band for stars in the range ~ 
04 to B7 (Figure ^ suggests a spectral type could be determined to ~ ± 2 sub-types). Note 
that He I 1.7002 /xm is detected in HD183143 (B7 la) even though the telluric correction 
appears worse for this star: compare the region just to the red wavelength side of the 1.7 
/im line. The dip in the spectrum is due to a relatively poor air-mass match between object 
and standard, but the He I 1.7002 /im line is visible in the uncorrected spectrum. 

Figure |^ shows a more detailed look at the ii— band spectra of the four O stars. The 
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figure suggests a tentative detection of He II (13-7) at 1.5719 /im and perhaps He II (12-7) 
at 1.6918 /im. The former is clearly in a region of higher signal-to-noise; the latter sits 
in a region with larger telluric features (the strongest of which are to the blue of the He I 
1.7002 /im line). The 1.5719 /im line appears in both of the O star spectra of Langon & 
Rocca-Volmerange (1992) as well, including their spectrum of HD 199579, one of the stars 
observed in our sample. Detection of He II is difficult since the 1.5719 /im line is near the 
Brackett 15-4 line (1.5701 /im); and hence, the He II line could depend on how well we have 
corrected the Brackett series lines in the A star which was used as a atmospheric standard 
(see §2). The He II lines may prove more useful when observed at higher spectral resolution. 



3.2. /' and J— band Spectra 

Figure ^ shows spectra of the OB stars near 1.1 and 1.3 /im. The primary lines which 
fall within our /' and J— band spectra are due to H and He I, as in the if —band spectra. 
However the H lines, Pa7 (1.0938 /im) and Pa/? (1.2818 /im) are blended with lines of He I 
and He II at our resolution (see Figure ^), so they are of limited use for determining spectral 
types. We note that He I 1.0830 /im is present in several of the stars (HD190429, P Cygni 
profile; HD 183143). We anticipate that the behavior of this triplet fine {2p^P°-2s'^S) will 
not be well correlated with spectral type, following the example of the analogous singlet 
He I feature at 2.0581 /im (HCR96). We have only five /' spectra in our initial sample, so 
we do not attempt to make any firm conclusions about this wavelength region. In any case, 
the /' and J— bands should be useful for comparison to model atmosphere predictions, 
especially when they are obtained for "free," i.e. in cross-dispersed mode. 



4. SUMMARY 

We have presented if —band spectra of a preliminary set of OB stars and identified 
absorption features that can be used to classify young, massive stars in obscured H II 
regions and/or along sight lines with large optical extinction due to interstellar dust. The 
absorption strength of the He I (4-3) 1.7002 /im line, in conjunction with the presence or 
absence of the if —band lines of the Brackett series, is well correlated with optical spectral 
types and can be used as a coarse spectral classification. The behavior of the infrared lines 
is analogous to well-studied lines in the optical wavelength region including basic effects 
due to excitation (temperature) and gravity (pressure). In particular, we see the effects of 
linear Stark broadening, resulting from changes in pressure in B dwarfs and supergiants, 
on the Brackett lines which will be useful as a luminosity indicator in the B stars. We 



have also presented spectra of the OB stars between 1.1 and 1.3 /im. Features in these 
wavelength regions may prove useful in cases where excess emission is still strong at H or 
for the purposes of testing atmosphere models of hot stars. 

Support for this work was provided by NASA through grant number HF 01067.01 - 
94A from the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. OSIRIS 
was built with support from NSF grants AST 90-16112 and AST 9218449. This work was 
also supported by NSF grant AST 93-14808. 
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Fig. 1. — A/AA ~ 570 if —band spectra of OB stars (The two latest B stars are from an 
independent data set and have A/AA ~ 525; see text). The spectra have been normahzed 
by a low order fit to the continuum. He I 1.7002 //m and H Brackett lines are the key to 
the present if —band classification of OB stars. The first Brackett line is BrlO. See text and 
Figure 0. The intensity scale is the same for each star. 



-12- 



4 — 



E 
o 



S 3h 
E 

o 



X 



(0 

> 

D" 



- « 

— o 

o o 

-To o T - 

^ o - 

- o ^ -L 

<> 



04 06 08 BO B2 B4 

Spectral Type 



B6 



Fig. 2. — He 1 1.7002 /xm equivalent width versus spectral type. O type stars and late B type 
stars are distinguished by the presence or absence of H Brackett absorption; see Figure |l]. 
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Fig. 3. — Detail of if —band spectra for the O stars (same spectra as plotted in Figure 0) 
showing probable He II absorption at 1.5719 /im. He II 1.6918 /xm absorption is marked 
although not clearly detected; see text. The intensity scale is the same for each star. 
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Fig. 4. — A/AA ~ 570 /' and J— band spectra of OB stars. The spectra have been normahzed 
by a low order fit to the continuum. Pa7 (1.0938 /im) and Pa/? (1.2818 fim) may have 
contributions due to He I and He II. The intensity scale is the same for each star. 



